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By Gaylord 0. Ellis and John D. Stanitz 


SUMMARY 

Two numerical examples are presented for two-dimensional, com- 
pressible, nonviscous flow in centrifugal compressors with backward- 
curved, logarithmic -spiral blades, the center lines of which (in the 
axial-radial plane) generate a right circular cone when rotated about 
the axis of the compressor. The two examples are for different blade 
angles and are compared with a third example for straight, radial 
blades. The solutions were obtained in a region of the compressors 
(including the impeller tip) that was considered to be unaffected by 
the impeller-inlet configuration and by the diffuser vanes. The 
numerical results are presented in tables of the stream-function 
distribution and in plots of the streamlines, constant Mach number . 
lines, and constant -pressure -ratio lines. In addition, a simplified 
analysis for logarithmic -spiral blades is presented that can be used 
to determine flow conditions within the impeller except near the 
impeller tip and the impeller inlet. 

For the same impeller-tip speed, number of blades, and compressor 
flow rate, it was concluded that: (1) The deceleration along the trail 
ing face of the blade (face opposed to the direction of rotation) is 
less for backward -curved blades than for straight blades; (2) the 
relative eddy, which forms on the driving face of the blade at low 
ratios of flow rate to impeller-tip speed, can be reduced or eliminated 
by backward -curved or forward -curved blades; (3) the static pressures 
are essentially equal at the impeller tip for straight and backward- 
curved blades. The slip factor, corrected for. the tangential component 
of the relative flow assuming perfect guiding of the fluid by the 
blades, was found to be independent of the blade angle for logarithmic - 
spiral blades. 
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INTRODUCTION 

Improved performance of centrifugal compressors should result from 
improved design methods based upon increased knowledge of the compli- 
cated flow processes in centrifugal compressors. A series of solutions 
have therefore been obtained at the NACA Lewis laboratory for two- 
dimensional, compressible, nonviscous flow through centrifugal com- 
pressors. Each solution shows the effect of varying a single design or 
operating variable from the fixed standard values common to the other 
solutions. The solutions are not intended to provide specific design 
information but are meant to clarify the flow phenomena in centrifugal 
compressors. 

A general method of analysis is developed in reference 1 for 
conical -flow compressors and turbines in which the center line of the 
passage generates a right circular cone about the axis of rotation. 

The radial -flow centrifugal compressor is a special case in which the 
cone angle is 180°. 

The general method of analysis is used in reference 2 to investi- 
gate the effect of compressor -flow rate, impeller-tip speed, and the 
number of blades upon the compressible flow within a certain region of 
radial- or conical -flow compressors with straight blades. The region 
investigated, which included the impeller tip, was assumed to be unaf- 
fected by the impeller-inlet configuration and by the diffuser vanes; 
that is, the impeller inlet and the diffuser vanes, if any, were con- 
sidered far enough removed from the region investigated not to appre- 
ciably affect the flow in that region. 

Solutions are extended herein to two impellers- with backward - 
curved, logarithmic -spiral blades. The solutions are obtained for flow 
regions of the compressor subject to the same restrictions outlined for 
the straight -blade solutions. The purpose of these solutions is to 
investigate the effect of logarithmic -spiral -blade curvature upon the 
flow conditions within the compressor, especially along the blade 
surfaces , and to obtain qualitative information on the effect of blade 
curvature in general. In addition, a simplified method of analysis is 
presented that permits the relatively rapid calculation of flow con- 
ditions within the impeller, except near the impeller tip and impeller 
inlet and which is used to draw some conclusions with regard to forward - 
curved logarithmic spiral blades. 
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. METHOD OF SOLUTION 

A general analysis is developed, in reference 1 for steady, two- 
dimensional, nonviscous, compressible flow in centrifugal compressors 
with arbitrary blade shapes and arbitrary variations in the passage 
height. The analysis is limited to radial- and conical -flow com- 
pressors in which the center line of the passage generates a right 
circular cone when rotated about the axis of the compressor (fig. l). 

The two-dimensional flow pattern is considered to lie on the surface 
of this cone. The final equations, developed from the equations of 
absolute irrotational motion, continuity, and energy, in reference 1 
are presented together with a brief discussion of the coordinate system, 
the assumptions and limitations of the general analysis, and the 
relaxation methods used in the numerical procedure. 

Coordinate system. - A developed view of the conic surface gen- 
erated by the passage center line (fig. l) is shown in figure 2. The 
dimensionless conic coordinates of a fluid particle on the .conic surface 
are R and 0. (All symbols are defined in the appendix. ) The conic - 
radius ratio R is defined by 



( 1 ) 


where r is the conic radius (distance along conic element from apex 
of cone) and the subscript T refers to the impeller tip. The 
coordinate system (R,0) rotates with the angular velocity CO of the 
impeller. The pas sage -height ratio H in the direction normal to the . 
conic surface (fig. 2) is a continuous function of the conic -radius 
ratio E. 

H = jjj » f(R) (2) 

where h is the passage height at any conic -radius ratio R. 

Assumptions and limitations. - This general analysis assumes that 
flow conditions are uniform across the passage normal to the conic sur- 
face; that is, the flow varies only along the conic surface. In order 
to satisfy this assumption, it is necessary that: (l) the gradient of 
h with respect to r be small; and (2) the cone angle a be suffi- 
ciently large. The allowable variation in a from 180° depends on 
the ratio h/r and on the desired accuracy. For the hypothetical 
'limiting case in which h/r approaches zero everywhere along the conic 
surface , the analysis is accurate for all values of a. 
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Velocity -ratio components. - The fluid particle on the developed 
conic surface in figure 2 has a relative tangential -velocity ratio U 
and a. radial- (along conic element) velocity ratio V. These velocity 
ratios are defined by 



o 


V « — (3b) 

c o 

where 

c local speed of sound 

u tangential component of velocity relative to impeller (positive in 
direction of rotation) 

v radial (along conic element) component of velocity 
Subscript: 

o absolute inlet stagnation condition 

Stream function. - A dimensionless stream function \J/ satisfies 
the continuity equation if defined as 


and 



(4a) 



(4b) 


where p is the weight density. 


The stream function \}/ is constant along blade surfaces . If 
is assigned the value 0 along the driving face of the blade (the blade 
surface in the direction of rotation), then the value of ^ along the 
trailing face of the blade (the blade surface opposed to the direction 
of rotation) is given by (reference 1) 
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where 


0 = 0 t"®d 


and where 


where 


cp 


V 

Wo 


( 6 ) 


(7) 


= BJ_r_h m 
T T T 


( 8 ) 


where 

O included passage angle between blade surfaces at constant radius 


<W 


cp compressor flow coefficient 
W i compressor flow rate 
a flow area (normal to conic surface) 
B number of passages ( or blades ) . 
Subscripts: . 


d driving face' of blade (blade surface in direction of rotation) 

t trailing face of blade (blade surface opposed to direction of 
'rotation) 

Differential equation for stream-function distribution. - The dif- 
ferential equation for the stream-function distribution in compressors 
with thin, logarithmic -spiral blades is given by equation (35) of 
reference 1 in terms of transformed coordinates | and tj 
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where the impeller-tip Mach number is defined as 


"r sin | 
o 


(9) 


( 10 ) 


and where P is the blade curvature angle between a radius (conic 
element) and the tangent to the blade surface (fig- 3). The angle P 
is a 'constant for a given logarithmic -spiral blade. For positive 
values of P, the blades are curved forward (in the direction of 
rotation); for negative values of P, the blades are curved backward; 
for P equal to 0, the blades are straight and lie along R. 

The transformed coordinates £ and q, which are related to 
R and 0 by (reference 1) 

£ = cp (log R + 0 tan P) , (11a) 

e 

q = cp (0 - tan P log R) (lib) 

© 

have been introduced because they result in straight, parallel blades 
in the transformed |q -plane. Such a transformation is desirable 
because it simplifies the solution of the differential equation by 
relaxation methods. Lines of constant £ and q in the physical 
R0-plane are shown for a typical logarithmic -spiral blade in figure 4. 
These transformed coordinates, given by equation ( 11 ) , are developed in 
reference 1 for the conditions that q = 0 along the driving face of 
the blade, q = \J/ = cpcL, along the trailing face, and 0=0 for 

q = 0 at the impeller tip R =1.0. 
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In eq.uati.on (9), the blade -height ratio H is assumed, for con- 
venience, to vary with the conic -radius ratio according to 

H =*E m (12) 

where m is an arbitrary exponent. (For m = 0 the blade height 
remains constant and for m = -1.0 the flow area remains constant.) 

In order to solve equation (9), it is necessary to know the 
density ratio, which is related to the impeller-tip Mach number and 
the relative velocity ratio Q by (equation (11), reference 1) 

1 • 

"f 1 + *ir (em t ) 2-q2 ] (13) 

where 7 is the ratio of specific heats and 

Q (l) 2 =u 2 +V 2 (14) 

and where the absolute whirl of the fluid ahead of the impeller is 
assumed to be zero. The velocity Q (multiplied by p/p ) is in 
turn given by equations (4), (11), (12), and (14) 0 

1 

(aF-ffiff 


Equations .( 9 ) , (13), and (15) provide three equations with three 
unknowns, p/p Q -, Q, and \J/. 

Numerical procedure. - The system of equations ( 9 ) , (13), and (15) 
is solved by relaxation methods to obtain the stream-function distri- 
bution within the compressor. From this distribution, the velocity 
components and other conditions can be determined using equations (4), 
and so forth. Detailed outlines for the numerical procedure are given 
in references 1 to 3. The procedure is briefly sketched herein. 

(l) Equations (9) and (15) are changed to finite -difference f orm. 

(reference 2). 
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(2) Values of \J/ are specified on the boundaries of the flow 
region. Because the flow region investigated by the solutions in this 
report doe3 not extend as far upstream as the impeller inlet, the 
boundary values of across the passage between blades at the upstream 
boundary are determined by the simplified analysis subsequently given. 
Other boundary values of \J/ are determined by methods similar to those 
in references 1 and 2. 

(3) Values of are estimated at equally spaced points of a grid 
system within the boundaries of the flow region. 

(4) The preceding estimated values of \U are adjusted (relaxed) 
by the relaxation process until they satisfy eq.uation (9) in finite - 
difference form. 

(5) The downstream boundary values of t in the vaneless diffuser 
are adjusted to satisfy the Joukowski condition for tangency of flow 
leaving the tip of the impeller blade. 

(6) After the Joukowski condition has been satisfied, the grid 
spacing is reduced near the impeller tip in order to obtain detailed 
knowledge of the flow characteristics in this region where conditions 
are rapidly changing. 


RELAXATION SOLUTIONS 
Numerical Examples 

Two numerical examples are presented for two-dimensional, com- 
pressible flow through impellers with backward-curved, logarithmic- 
spiral blades with the blade angles 8 (fig* 3) equal to tan"l(-’6.5) 

and tan '*'(-1.0). A third example, for straight blades 8 = tan ^(O), 
is also given. Logarithmic -spiral blades were selected because the 
known transformation of coordinates given by equation (11) for this 
particular blade shape considerably reduced the labor that would 
have been involved in the relaxation solution of an arbitrary blade 
shape. It is considered that the results of these solutions can be 
applied in a qualitative manner to blades with arbitrary curvature., 
The other design and operating parameters, which are the same for all 
three examples, are as follows: 
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Flow coefficient, cp 0.5 

Impeller-tip Mach number, M^, • 1.5 

Exponent for variation in H with R, m -1.0 

Ratio of specific heats, y 1.4 

Included passage angle between blade surfaces 

at a constant radius, O , deg 18 


(For infinitely thin blades, a is independent of R so that. 

Qj, = C = constant. ) 

The solutions were . obtained in a region of the compressors 
(including the impeller tip) that was considered to be unaffected by 
the impeller-inlet configuration and by the diffuser vanes; that is, 
the diffuser vanes (if any) and the impeller inlet must be far enough 
removed from the region investigated not to affect appreciably the 
flow in that region. Each solution applies , within the limitations 
imposed by the assumption of two-dimensional flow, to conical -flow com 
pressors (and turbines) with various cone angles a (including the 
radial -flow compressor, a = 180°) but with the same included passage 
angle a (reference l). 

Results 

The numerical results are presented in plots of the streamlines , 
constant Mach number lines, and constant -pressure -ratio lines in fig- 
ures 5 to 7. The value of obtained by relaxation methods at each 
of the grid points in the flow field is given for the three numerical 
examples in tables I to III. The velocities, pressures, and so 
forth, can be obtained from these values of ii with equations (4), 
(13), and (14) by methods outlined in reference 1. The information 
obtained from these tables could be used to develop correlation 
equations (similar to those in reference 2) for predicting flow con- 
ditions in the impeller at other, design and operating conditions, 
although past attempts to develop such correlation equations for 
logarithmic -spired blades have failed. 

Streeml ines . - The streamline configurations (relative to the 
impeller) for the three examples are shown, in figure 5. In order to 
assist in the proper geometrical orientation of the curved-blade 
examples, a dashed outline of the straight -blade impeller is superim- 
posed on the flow field of the curved -blade solutions in figures 5 to 
7. The streamlines are designated in such a manner (\|//^) that the 
value of a streamline indicates the percentage of the total flow (in 
the passage) between .the streamline and the driving face of the blade. 
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For a given density ratio and blade height, the streamline spacing is 
indicative of the velocities relative to the impeller, with close 
spacing indicating high velocities and wide spacing indicating low 
velocities. 

The streamlines for the straight -blade example are given in fig- 
ure 5(a). For the design and operating conditions of this example, an 
eddy has formed on the driving face of the blade. This eddy indicates 
negative velocities on and near the blade face and results from the 
condition of absolute irrotational motion. The eddy is attached to 
the blade and rotates with an angular velocity equal and opposite to 
the rotational velocity of the impeller. If the flow rate through the 
impeller were 0, this eddy would occupy the entire passage. In actual 
practice this eddy is probably unstable and should be, eliminated by 
proper changes in the design and operating conditions of the compressor. 
From a comparison of figures 5(a), 5(b), and 5(c), it appears that the 
eddy can be reduced or eliminated by backward-curved logarithmic -spiral 
blades (for the same impeller-tip speed, number of blades, and compres- 
sor flow rate). 

Mach number lines. - Lines of constant Mach number relative to 
the impeller' are shown for the three examples in figure 6 . The gen- 
eral characteristics of these plots are similar. The velocities (as 
indicated by the Mach number lines) along the driving face of the 
blade are low; the velocities along the trailing face are high; and 
the velocities become equal on both faces at the blade tip (as required 
by the Joukowski condition). The maximum relative Mach number occurs 
on the trailing face of the blade at a radius ratio well within the 
impeller and the flow decelerates along the face of the blade from 
this point to the blade tip. This deceleration is seen to decrease, 
however, as the blades become more backwardly curved. 

If boundary -layer effects are neglected, the nonuniformity of 
the velocity at the impeller tip becomes negligible at a radius ratio 
of approximately 1.15. The average relative Mach number at the impeller 
tip is low because of the high impeller-tip Mach number, which results 
in high fluid densities and therefore low velocities. 

Pressure -ratio lines. - Lines of constant static -pressure ratio 
(local pressure divided by the absolute inlet stagnation pressure) are 
shown for the three examples in figure 7. The general characteristics^ 
of these plots are the same. At a given radius ratio, the pressure 
ratio is higher on the driving face of the blade than on the trailing 
face, except at the blade tip where the pressure ratios are equal. 

This difference in pressure ratio accounts for the impeller torque. 
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It is interesting to note that the static -pressure ratios are essen- 
tially equal at the impeller tip in all three examples, although the 
total pressure must he progressively less as the blades become more 
backwardly curved. At the impeller tip, the absolute kinetic energy 
of the fluid (which is measured by the difference in total and static • 
pressures) therefore becomes progressively less as the blade becomes 
more backwardly curved. The diffuser problem should therefore be less 
critical for backward -curved blades than for straight blades. 

Velocity ratios along blade surfaces. - Of special interest, 
because of boundary -layer considerations, are the relative velocity- 
ratio distributions along the trailing and driving faces of the impeller 
blades, and Q^. These velocities are plotted in figure 8 as a 

function of the radius ratio E. Negative values of Q along the driving 
face of the blade indicate the presence of a relative eddy on the blade 
surface. The deceleration is seen to be less rapid on the trailing 
face of the blade near the tip when the blade is curved farther back- 
ward. The velocities along the trailing face of the blade away from 
the impeller tip (E = 1.0 at tip) are about equal for all three 
examples. This equality results because, although the difference in 
velocities on the two blade surfaces becomes less as the blades are 
curved farther backward, the mean relative velocity must increase to 
accomplish the required through flow, and as a result the net change in 
velocity on the trailing face with changes in. P is negligible (for 
the range of P investigated in the three numerical examples ) . 

A measure of the tendency for boundary -layer separation is given 
by the term 


8 dQ ■ 

Q ds 

where 6 is the boundary -layer thickness and s is distance along the 
blade surface. In the case of laminar flow, for example, this term must 
not become more negative than a certain value (the magnitude of which 
is not discussed in this report) or separation will occur. In dimen- 
sionless form this term becomes 
✓ 

s / r T ag 

Q as 

where S is distance along the blade surface in units of the tip radius. 
For a given value of &/r T , the more negative the value of ^ the 
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more likely will boundary -layer separation occur. In figure 9, 

i. $ 0 - is shown as a function of R at the critical region on the trail - 
Q dS 

ing face of the blade near the impeller tip for the three numerical 
examples. It is evident that backward-curved logarithm! c -spiral 
blades are more desirable than straight blades with regard to velocity 
gradient on the trailing face of the blade (at least for the same 
impeller-tip speed, number of blades, and compressor flow rate). If 
it is desired, to compare the velocity gradients for straight and 
logarithmic spiral blades at the samfe work input, that is, at the 
same value of pM^, < (where the slip factor p is defined as the ratio 

of the average absolute tangential velocity of the fluid leaving the 
impeller tip to the tip speed of the impeller), then 

(PJV 2 ) =(^n 2 ) (16) 

1 straight log -spiral 

For the numerical 'example with 3 equal to tan~^(-1.0), the 
"equivalent" tip Mach number for a straight blade is, from equation (16), 


1 S xi - 5 =i - 3868 

where the numerical values for the slip factor are given in the next 
section. (Note that for straight blades the slip factor is independent 
of impeller-tip speed, reference 2. ) With the use of the correlation 
equations of reference 2, the velocity gradients were determined for 
straight blades at the equivalent tip Mach number of 1.3868 and the 
result is given by the dashed curve in figure 9. For straight blades 
the velocity gradients were changed very little by the reduced 
equivalent tip Mach number so that the velocity gradients for 
logarithmic -spiral blades are definitely more desirable than those 
for straight blades at the same work input. 

Slip factor. - A method of evaluating the impeller slip factor 
p is given in appendix E of reference 1. The following values of p 
were obtained for the three* examples of this report: 


Blade angle 
P 

Slip factor 

tan“^(0) 

0.8986 

tan”^( -0.5) 

.8343 

tan 1 ( -1.0) 

.7681 
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As would be expected, the slip factor is less for backward -curved 
blades than for straight blades because for backward -curved blades 
the tangential component of the relative velocity is opposed to the 
direction of rotation. If the slip factor is corrected by adding 
this tangential component of the relative flow, assuming the fluid is - 
perfectly guided, the corrected slip factors become 


Blade angle 

Slip factor 

0 

^(corrected) 

tan”^(0) 

0.8986 

tan -1 ( -0.5) 

.9013 

tan -1 ( -1.0) 

.9036 


These corrected values of p are essentially equal and it is therefore 
concluded that the slip velocity, which is defined as the average 
absolute tangential velocity at the impeller tip minus the absolute 
tangential velocity with perfect guiding by the blades , that is , 
assuming an infinite number of blades, is the same for various negative 
values of blade angle 0. This conclusion is not in agreement with 
the well-known Stodola correction in which the slip velocity varies 
directly with cos 0 (reference 4, p. 1259). 


SIMPLIFIED ANALYSIS 

The relaxation solutions presented herein require a large number 
of calculations . It would therefore be advantageous to have a quicker 
(although less accurate) means of estimating the flow conditions in 
impellers with logarithmic -spiral blades. In this section, a simplified 
method of analysis is developed whereby flow conditions within the 
impeller, except near the tip, can be estimated for radial- and conical - 
flow impellers with logarithmic -spiral blades. The simplified analysis 
is based on the assumption that for high-solidity blades, such as those ' 
in centrifugal compressors , the component of the relative velocity 
normal to the blade is 0 at all radii within the impeller; that is, 
the fluid is perfectly guided by the blades. 

Velocity distribution. - In terms of U, V, and M^, the condition 
of irrotational motion given by equation (5) of reference 1 becomes 
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ou u . Su ■ 1 dv 

' '“'i i’5'55 

But, from the assumption of this simplified analysis, 

U = Q sin 3 
y = ft cos 3 

so that equation (17) becomes 

= sin 3 § + sin 3 if 


-2^/ 


R 


cos 3 oft 
r 5e 


In terns of the transformed coordinates £ and rj, given by 
equation (11) and shown as functions of R and 0 in figure 4, 
equation (18) becomes 

(i - q tan 3\ ^ „ cp dft 

-2Hp 2 - « * ln p - 351 dwi 

\ cp sec 3 

which integrates to give 

^ exp ( | ~ T ^ an A - constant exp ( T ^ an . 


sin 3 


\cp sec 3 


yep sec 3y 


The constant of integration is determined from the condition that 
Q = ft d when i) = 0 (fig. 4); equation (19) then becomes 


(17) 


(18) 


(19) 


ft = ft d exp 


f r\ tan 3 
lep sec 2 3/ 


2Mm 

+ si5T elp 


'—hr ) 811,11 

icp sec 3/ \cp sec 3 j 


( 20 ) 


Equation (20) gives the variation in ft with t) along a line 
of constant £. (See fig. 4.) The determination of ft d from con- 
tinuity considerations is given in the following section. 


If 3~ >0, then from equation (11) £ = cp log e R • and t| = cp0; 

therefore, after- the indeterminates are evaluated, equation (20) 
becomes 


ft = ft d + SRM^ 


(20a) 
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which agrees with equation (40) of reference 1 for straight blades 
(3=0). 

The value of £ in equation (20) can be expressed in terms of the 
radius ratio R^ at which i cuts the driving face of the blade 

(q = 0). From equation (11) with R = R^ and q = 0, 

| = cp(log e R d + 0 tan 3) 

0 = <p(0 - tan 3 l°S e R^) 
so that, after eliminating 0, 

R = exp [ ^-5 — j (21 

a \cp sec 3/ 

which can be substituted into equation (20). 


Stream function. - The total differential for the 'stream function 


is given by 


^ ■ sl de + Ir “ 


which, from, equation (4), becomes 

d\J/=- e -VHRd0 - — UH dR 


p o 


= 


(• 


HRQ ( cos 3 d0 - sin 3 


dR 
E > 


( 22 ) 


But, from equation (11), 

dq = cp(d0 - tan 3 
so that equation (22) becomes 


= cp(d< 


dR 
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cL^ = -2- HBQ cos 6 (23) 

p o 9 

where p/p Q is given by equation (13). 

In order to determine Q d in equation (20), equation (23) is 
numerically integrated across the passage to give 

A * +t * cpo T 

i-BE«dn (24) 

Jo ° 

which equation is satisfied when the proper value of Q d is con- 
tained in the expression for Q (equation (20)). After the value of 
Q d is determined, the variation in \{/ with q across the passage can 
be determined by the numerical solution of equation (23). 

Comparison with relaxation solutions. - Comparisons between the 
results of the simplified analysis and the relaxation solutions are 
given in figures 10 to 12. In figure 10 the velocity ratio along the 
driving face of the blade Q d is plotted as a function of the radius 

ratio along the driving face of the blade B^ for the three values of 

6: tan’^O), tan _1 (-0.5), and tan" 1 (-1.0). The results of the 

simplified analysis are shown by the heavy curves; the relaxation 
solution is indicated by the plotted points. The agreement is con- 
sidered to be good up to a conic radius ratio of 0.800 or 0.825. (This 
agreement exists for all values of B< 0.800 provided the impeller inlet 
can be neglected. ) 

In figure 11 the velocity ratio along the trailing face of the blade 
is plotted as a function of the radius ratio along the trailing face 

of the blade B for the three values of 0 used. in the numerical 
examples. (Substantial agreement between the simplified analysis and 
the relaxation solution exists for all values of B less than those 
shown on the plot . ) The value of B^ at which agreement fails , depends 

on the value of 0. For 0 equal to 0, the agreement is satisfactory 
up to B t equal to 0.800 or 0.825; for 0 equal to tan -1 (-0.5) the 
agreement is satisfactory up to B^ equal to approximately 0.9; for 
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3 equal to tan"^( -1.0) the agreement is satisfactory up to R equal 

"G 

to approximately 0.95. From equation (11), it can he shown that these 
values of R^ correspond to the values of R^ (for the same values of 

| ) at which the agreement is satisfactory for along the driving 

face of the blade (R^ = 0.800 or 0.825, fig. 10). The ratio of 
R^/R^ for a given value of £ is developed later and presented in 
equation (28). For 0^, equal to 18° and for the range of 3 investi- 
gated, the agreement between the relaxation solution and the simplified 
analysis therefore satisfactory for values of i corresponding to 
values of R^ up to 0.800 or 0.825. 

In figure 12 the velocity ratio Q is plotted across the passage 

along lines of constant | for three values of. 3. The value of i 

for each value of 3 was selected to have a value of R, between 

d 

0.800 and 0.825. The agreement between the relaxation solution and 
the simplified analysis is satisfactory. 


Particular set of solutions. - In general, equation (24) must be 
integrated by numerical methods. If, however, the fluid is incom- 
pressible (p/p =1) and if the flow area is constant (HR = l), 
equation (24) combined with equation (20) integrates across the passage 
from T) and = 0 on the driving face of one blade to rj and \| / ='!' 
on the trailing face of the next blade to give 



which, combined with equations (5) and (21), becomes 



sin 3 

sec 3 / 


O £2 — 

™T. 


^d 

sin^ 3 


cosh 


'0_ sin 3 

_T 

sec 3 



( 25 ) 


Equation (25) gives the variation in C^/M^, with 

cp/M^. For incompressible flow, the speed of sound contained in the 


R d’ °T y and 
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denominators of Q, Hp, and cp has no. physical significance and may 

have any arbitrary value because these quantities appear as ratios in 
equation (25). The quantity is the ratio of the velocity along 

the driving face of the blade to the tip speed of the impeller and 
cp/Hj, is proportional to the ratio of the flow rate through the impeller 

to the tip speed of the impeller. It will be noted in equation (25) 
that ft^Mp decreases linearly with increasing radius ratio R^. 

The ratio Q^/Hp along the trailing face of the blade can be 

determined from equation (20) with q =c PO^, and with R^ given by 
equation (21). • 


Hp 


2R d /a T sin 0\ Q d /0 T sin 0 

sin~0 Sinh \ sec V) + eXp \TZT 


(26) 


where Q^/Hp is gl ven equation (25).. The value of R d contained 

in equations (25) and (26) is related to the value of R ^ , at which 
occurs, as follows: 

From equation (11) eliminating 6 


„ ft - r\ tan 0' 

R = exp — 

\.cp sec 0 

From equation (21) and because R = R^. when q =cpa^,, equation (27) 
becomes 


(27) 


R t = R d exp 


(~ 0 T sin 0V 
1 sec 0 / 


(28) 


Equation (28) relates R^ and R d the ratio of which is seen to be 

constant for a given impeller with logarithmic -spiral blades. It will 
be noted from equation (26) and (28) that ft t Mp increases linearly 
with increasing values of R^. 

In order to evaluate in a qualitative manner the effect of 0 
upon the values of Q d and the variations in Q d /Hp and 

with 0 for various values of R are plotted in figure 13 from 
equations (25) and (26) for 0^, = 0.31416 (18°) and cp/M,^ = 0.250. It 

will be noted that 0 has a pronounced effect upon both Q^/Hp a 11 *! 
for large absolute values of 0. In particular the value of 
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increases at an increasingly greater rate as the absolute value 

increases from 0. Because the minimum value of Q^/Hp occurs 

in the immediate neighborhood of 3=0, the tendency for an eddy to 
form on the driving blade surface is greatest for straight blades. 
Therefore, at the same impeller tip speed, number of blades, and com- 
pressor flow rate, forward-curved blades (3>0), as well as backward- 
curved blades (3<0), have less tendency than straight blades (3 = 0) 
for the formation of undesirable relative eddy flows on the driving 
face. This fact might be unexpected because for forward -curved blades 
the loading is higher and therefore the velocity on the driving face 
might be' expected to be lower. In order to accomplish the required 
through flow, however, the relative velocities must increase as the 
blades are curved forward. This fact explains the very pronounced 
effect of large absolute values of 3 upon Q /M^, and q /n^. 


of 3 


In figure 13 the variation (with 3) in the relative velocity 
ratios at the driving and trailing faces of the blades and Q t 

at a radius ratio of 0.75 (for example) agrees with the variation shown 
in figure 8. In particular, in both figures the value of is some- 
what lower for 3 = tan" -*-(-0.5) than for tan" ■*-(-!. 0) or tan" ^ ( 0 ) . 


SUMMARY OF RESULTS AND CONCLUSIONS 

Two numerical examples are presented for two-dimensional, com- 
pressible, nonviscous flow in centrifugal compressors with backward- 
curved, logarithmic -spiral blades, the center lines of which (in the 
axial -radial plane) generate a right circular cone when rotated about 
the axis of the compressor. The two examples for spiral blades are 
for different blade angles and are compared with a third example for 
straight, radial blades. All other design and operating parameters are 
the same for the three examples. The solutions were obtained in a 
region of the compressors (including the impeller tip) that was con- 
sidered to be unaffected by the inlet configuration of the impeller 
and by the diffuser vanes. (That is, the impeller inlet and the 
diffuser vanes, if any, were far enough removed from the region 
investigated not to appreciably affect the flow in that region. ) Each 
solution applies to radial- and conical -flow compressors (and turbines 
with. the rotation and flow direction reversed) with various cone angles 
but with the same included passage angle between blades. The numerical 
results are presented in tables of the stream-function distribution and 
in plots of the streamlines, constant Mach number lines, and constant - 
pressure -ratio lines. 
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In addition, a simplified analysis for logaritnmic-spiral blades 
is presented that can be used to determine flow conditions within the 
impeller except near the impeller tip and impeller inlet. For the 
range of blade angle 3 investigated and for an included passage angle 
between blades of 18°, the agreement between the relaxation solution 
and the simplified analysis was satisfactory for values of transformed 
coordinate | corresponding to values of conic -radius ratio along the 
driving face Q of blade up to 0.800 or 0.825. For a backward -curved 
blade of -45 and an included passage angle between blades of 18° 
the corresponding value of conic -radius ratio along' the trailing face 
of blade at which agreement was good is as high as 0.95. 

The principal conclusions of the work presented herein are: 

1. The maximum relative Mach number occurs on the trailing face 
of the blade at a radius ratio well within the impeller and the flow 
decelerates along the face of the blade from this point to the blade 
tip. For the same impeller-tip speed, and so forth, this deceleration 
becomes less as the blades become more backwardly curved. Because of 
this lower rate of deceleration, backward -curved blades are more 
desirable than straight blades as far as boundary -layer separation is 
concerned. 

2. The relative eddy, which appears on the driving face of the 
blade at low ratios of flow rate to impeller-tip speed, can be reduced 
or eliminated (at the same impeller-tip speed, and so forth) by 
backward -curved logarithmic -spiral blades. (Also see conclusion 5.) 

3. The static -pressure ratios are essentially equal at the impeller 
tip for straight and backward -curved, logarithmic -spiral blades at the 
same impeller-tip speed. This fact indicates that at the impeller out- 
let the absolute kinetic energy of the fluid becomes progressively less 
as the blades become more backwardly curved. The diffuser problem 
should therefore be less critical for backward -curved blades. 

4. The slip factor (defined as the ratio of the average absolute 
tangential velocity of the fluid leaving the impeller tip to the tip 
speed of the impeller) decreases as the blades are curved farther back- 
ward. If the slip factor is corrected by accounting for the tangential 
component of the relative flow, however, assinning perfect guiding of 
the fluid by the blades, this corrected slip factor is essentially 
independent of the blade angle. 

5. From the simplified analysis, it is concluded that (at the 
same impeller-tip speed, and so forth) forward -curved blades (3 > 0), 
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as well as_ backward- curved blades (0<O), have less tendency than 
straight blades (0 = 0) for the formation of undesirable relative 
eddy flows on the driving face of the blade. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, August 7, 1950. 
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APPENDIX - SYMBOLS 


a 

B 

c 


exp 


H 

h 


The following symbols are used in this report: 
flow area normal to conic surface 
number of passages (or blades) 
local speed of sound 
exponential, jexp(x) = e x J 
passage -height ratio, h/h^ 

passage height normal to conic surface 

0)r T sin | 

impeller-tip Mach number, 

c 

o 


m passage -height exponent 

p absolute static pressure 

I 2 2 ~ 

Q relative velocity ratio , V U +V 

R , conic-radius ratio , r/r^ 

r conic radius (distance along conic element from apex of cone) 

S dimensionless distance along blade surface, s/r^ 

s distance along blade surface 

U relative tangential -velocity ratio, u/c q 

u relative tangential velocity 

V radial -velocity ratio, v/c q 

v radial velocity along conic element 

W total compressor flow rate 
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a cone angle (fig. l) 

P blade angle (fig. 3) 

c 

y ratio of specific heats 
6 boundary -layer thickness 
i) transformed coordinate, equation (lib) 
d angle, radians unless • otherwise specified 
H impeller slip factor 

% transformed coordinate, equation (lla) 

p weight density of fluid 

C included passage angle between blade surfaces at constant radius, 
e t" e d 


cp 


compressor flow coefficient , 


Wo 


\J/ compressible stream function 
Gt) impeller angular velocity 
Subscripts: 

d driving face of blade (blade surface in direction of rotation) 
o absolute inlet stagnation condition 

T impeller tip 

t trailing face of blade (blade surface opposed to direction of 
rotation) 
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TABLE I - DISTRIBUTION OF STREAM FUNCTION Y IN £ i] -PLANE FOR p = tan" 1 (0) 

[Flow coefficient cp, 0.5; impeller-tip Mach number , 1.5; exponent . m (for 
variation in H with R), -1.0; included passage angle between blades O, 18* 
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Figure 4. - Lines of constant j; and T) (transformed coordinates) for 
logarithmic-spiral blades in physical R9-plane. 
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Figure 5- - Relative streamlines for compressible flow- 
through centrifugal impellers with various values -of. 
blade angle P. Streamline designation i f indicates 

percentage of total flow (through passage) between 
streamline and driving face of blade (right side of , 
passage,). Flow coefficient cp, 0.5; impeller-tip Mach 
number Mj, 1.5; exponent m (for variation in H with 
R), -1.0; included passage angle between blades o, 18°; 
compressible flow (7 = 1.4). 
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Relative 

streamlines 


ig. 5(a) 


Angle, 0, deg 

(t>) Backward-curved blades, p * tan -1 (-0.5). 

Figure 5* - Continued. Relative streamlines for compressible 
flow through centrifugal impellers with various values of the 
•blade angle p . Streamline designation fty t indicates percentage 
of total flow (through passage) between streamline and driving face 
of blade (right side of passage). Flow coefficient Cp, 0.5; impeller- 


tip Mach ndmber 


1.5; exponent m (for variation in 


-1.0; included passage angle between blades 0. l8°: compressible flow 

( 7 = 1 .*). 
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(c) Backward- curved blades, p = tan" 1 (-1.0). 


Figure 5. - Concluded. Relative streamlines for compressible flow through centrifugal 
impellers with various values of the blade angle p . ' Streamline designation 
indicates- percentage of total flow (through passage) between streamline and driving 
face of blade (right side of passage). Flow coefficient <£,'0.5; impeller-tip Mach 
number 1.5; exponent m (for variation in H with R), -1.0; included passage 

angle between blades G , 18 °; compressible flow (7 » 1.4). 
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Figure 6. - Concluded, Lines of constant Mach. number relative to impeller. Flow 

coefficient cp, 0.5; impeller-tip Mach number M^, 1.5; exponent m (for variation in 
H with R ) t -1.0; included passage angle between blades O, 18°; compressible flow 
(7 ° 1.^). 
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Angle, 9, deg 

(b) Backward-curved blades, f3 * tan** 1 (-0.5). 

Figure 7* - Continued. Lines of constant static -pressure ratio. 
Flow coefficient cp, 0.5;. impeller- tip Mach number Mr, 1.5; 
exponent m (for variation in H with R), -1.0; included 
passage angle between blades C, 18°; compressible flow 

(y = i.*0. 
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Figure 7. - Concluded. Lines of constant static-pressure ratio. Flow coefficient cp, 0.5; 
impeller-tip Mach number Hy, 1,5; exponent m (for variation in H with R), -1.0; 
included passage angle between blades O, l8°; compressible flow (7 = 1.4). * 
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Figure 8. - Variation in relative ratio along driving and trailing faces of 

logarithmic -spiral blades for three values of blade angle f3 . Flow coefficient 
cp, 0.5; impeller- tip Mach number Mt, 1.5; exponent m (for variation in H 
with R), -1.0; included passage angle between blades O , 18°; compressible 
flow (7 = 1.4). 
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Conic-radius ratio along driving face t of blade, 

Figure 10. - Comparison of variation in relative velocity ratio along driving face of 
blade with conic-radius ratio along driving face of blade as obtained by relaxation 
methods and by the simplified analysis for three values of blade angle g. Flow 
coefficient cp, 0.5; impeller-tip Mach number M^, 1.5; exponent m (for variation 
in H with R), -1.0; included passage angle between blades 0 , 18°; compressible 
flow (7 = 1.4). ‘ ' . - 
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Figure'll. - Comparison of variation in relative velocity ratio along trailing face of 
blade with conic-radius ratio along trailing face of blade as obtained by relaxation 
methods and by the simplified analysis of three values of blade angle p . Flow 
coefficient cp, 0.5; impeller-tip Mach number Mij>, 1.5; exponent m (for variation 
in H with Ft), -1.0; included passage angle between blades a, 18°; compressible 
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Figure 12. - Comparison of variation in relative velocity ratio with ratio 
of distance across passage (along lines of constant |) as obtained by 
relaxation methods and by the simplified analysis for three values of 
blade angle j3 . Flow coefficient cp f 0.5; impeller-tip Mach number 
Mrp f 1.5; exponent m (for variation in H with R), -1.0; included passage 
angle between blades a, 18 °; compressible flow (7 = 1.4). 
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Figure 13. - Variation in velocity ratio at driving and trailing faces of 
blade with blade angle for several values of conic-radius ratio. Ratio 
cp/Mt, 0.250; exponent m (for variation in H with R) , -1.0; included 
angle between blades across the passage Op, 18°; incompressible flow. 
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